TECHNICAL  REPORT  1828 
September  2000 


The  Extended  COTS  LPI 
Communication  System 

Phase  3:  Performance  Analysis 

j.  c.  Allen 
M.  Reuter 
G.  M.  Dillard 
G. Johnson 
J.  Zeidler 


Approved  for  public  release; 
distributiori  is  unlimited, 


SSC  San  Diego 


20001030  138 


TECHNICAL  REPORT  1828 
September  2000 

The  Extended  COTS  LPI 
Communication  System 

Phase  3:  Performance  Analysis 


J.  c.  Allen 
M.  Reuter 
G.  M.  Dillard 
G.  Johnson 
J.  Zeidler 


Approved  for  public  release; 
distribution  is  unlimited. 


SSC  San  Diego 
San  Diego,  CA  92152-5001 


Irrio  II'SrerBOllBD  4 


SSC  SAN  DIEGO 

San  Diego,  California  92152-5001 


Ernest  L.  Valdes,  CAPT,  USN  R-  C.  Kolb 

Commanding  Officer  Executive  Director 


ADMINISTRATIVE  INFORMATION 

The  work  described  in  this  report  was  performed  for  the  Office  of  Naval  Research  by  the  Applied 
Research  &  Technology  Branch  (D363),  SSC  San  Diego. 


Released  by 

G.  W.  Anderson,  Head 

Applied  Research  &  Technology 

Branch 


Under  authority  of 
R.  H.  Moore,  Head 
Environmental  Science: 
Division 


SB 


Abstract 


Phase  3  of  the  COTS  LPI  effort  examines  two  LPI  systems:  LPI-1  approximates  an 
incoherent  M-ary  orthogonal  system;  LPI-2,  developed  by  Johnson  [19],  is  a  variant 
of  the  cychc  code-shift  keying  (CCSK)  waveform.  Each  LPI  system  is  assessed  using 
the  performance  vector  [Pg  Ei/No  R/W  SNR]^.  The  SNR  is  included  with  the  three 
standard  COMM  variables  as  a  measure  of  covertness.  Slices  of  this  function  are 
used  to  rank  LPI-1  and  LPI-2.  For  additive  white  Gaussian  noise,  the  fundamental 
trade-oflf  that  emerges  for  fixed  bandwidth  W  is  between  the  SNR  and  the  bit  rate  R. 
A  more  covert  system  requires  longer  symbol  vectors.  These  longer  symbol  vectors 
lower  the  bit  rate.  The  LPI  performance  sinrface  permits  system  engineers  to  graph¬ 
ically  assess  the  trade-oflfs.  For  the  multipaths  of  the  Extended  Littoral  Battlefield, 
simulations  show  that  multipath  may  either  have  no  effect,  improve  performance, 
or,  degrade  performance.  The  basic  observation  is  the  longer  symbol  vectors  can 
mitigate  the  multipath  effects.  Analysis  of  this  multipath  effect,  increasing  the  bit 
rate,  and  migrating  these  results  to  an  LPI  performance  surface  are  the  basic  Phase 
4  recommendations. 
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The  COTS  LPI  Communication 
Systems 

1.1  Migrating  COTS  to  LPI 

The  goal  of  this  project  is  to  migrate  commercial  off-the-shelf  (COTS)  technologies 
into  military  low-probability  of  intercept  (LPI)  applications.  This  task  was  under¬ 
taken  by  SPAWAR  Systems  Center,  San  Diego  because  commercial  wireless  companies 
cannot  profitably  invest  research  and  development  dollars  for  this  mihtary  market. 
Phase  1  of  this  effort  started  by  evaluating  COTS  cellular  and  COTS  Wireless  Lo¬ 
cal  Area  Networks  for  the  Marine  Corps  Littoral  Warfare  Applications  [14].  This 
evaluation  showed  that  several  COTS  systems  can  be  effective  in  many  mihtary  ap¬ 
plications  with  only  minor  modifications.  However,  these  systems  did  not  provide 
LPI  communications  for  some  httoral  operations. 

Phase  2  showed  how  to  overcome  these  hmitations  [19]:  “components  developed 
for  use  by  the  COTS  PCS  systems  could  be  used  in  a  communications  apphque  to 
provide  LPI  capability  far  beyond  that  of  the  original  COTS  system.”  The  successful 
construction  and  deployment  of  a  prototype  LPI  system  demanded  that  this  LPI 
capability  be  exploited. 

Phase  3  undertakes  this  task  to  (1)  model  the  LPI  receiver,  (2)  analyze  the  ba¬ 
sic  LPI  system  performance,  (3)  assess  LPI  system  performance  in  credible  radio¬ 
frequency  (RF)  environments.  This  report  documents  this  Phase  3  effort  and  lays 
the  foundation  for  the  Phase  4  effort. 

Section  2  sets  up  a  framework  to  analyze  the  LPI  systems  at  baseband.  Although 
the  LPI  waveform  is  not  new,  its  implementation  is  novel  and  several  LPI  variants 
readily  drop  out  of  this  framework.  We  analyze  two  LPI  systems  that  delineate 
performance  extrema:  LPI-1  is  an  approximation  to  an  M-ary  orthogonal  system, 
LPI-2  is  a  variation  of  cyclic  code-shift  keying  (CCSK).  Basic  to  any  communication 
system  is  its  probability  of  error.  This  section  obtains  the  probability  of  symbol  error 
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Pe  in  additive  white  Gaussian  noise  (AWGN)  and  makes  selected  trade-offs  with 
Eb/No  and  the  spectral  efficiency  R/W  [34].  However,  the  need  for  covertness  places 
additional  constraints  on  the  system.  A  weU-known  and  fundamental  variable  for 
covert  RF  operations  is  the  signal-to-noise  ratio  (SNR).  When  trades-offs  are  made 
between  the  probability  of  error  Pe,  the  spectral  efficiency  R/W,  and  the  SNR,  what 
emerges  is  a  two-dimensional  surface  in  Pe  x  R/W  x  SNR  space.  We  expect  that  RF 
mission  planners  will  exploit  this  novel  LPI  performance  surface. 

Communication  systems  do  not  operate  in  white  noise  but  must  successfully  over¬ 
come  their  RF  environments.  For  example,  anecdotes  from  the  Persian  Gulf  describe 
how  communication  platforms  could  be  seen  but  RF  contact  could  not  be  made. 
Some  LPI  operations  are  the  high-risk,  covert  forward  deployments.  It  is  uncon¬ 
scionable  to  put  our  forces  in  such  high-risk  deployments  without  proper  RF  mission 
planning.  More  generally,  we  cannot  simply  build  an  LPI  system  and  hand  it  off 
without  wrapping  the  system  in  an  operational  context.  The  trade-offs  of  bit  rate, 
error,  and  covertness  in  the  presence  of  multipath  are  non-trivial.  For  this  reason. 
Section  3  analyzes  the  Extended  Littoral  Battlefield  (ELB).  The  short-term  objective 
is  to  have  a  suite  of  credible  multipath  scenarios  to  assess  the  LPI  systems.  The 
long-term  objective  is  to  understand  RF  propagation  effects  in  the  ELB  and  then 
bring  this  knowledge  into  RF  mission  planning  and  LPI  deployments. 

Section  4  uses  the  preceding  results  to  assess  LPI  performance  in  the  ELB  multi- 
path.  ELB  multipath  can  either  degrade  or  improve  LPI  performance.  For  example, 
selecting  the  longer  partial  sequences  for  LPI-2  permitted  it  to  overcome  the  multi- 
path  degradations  in  the  HMMWV  scenario.  In  fact,  the  longer  sequences  eventually 
permitted  LPI-2  to  beat  the  AWGN  performance.  In  contrast,  it  was  the  shorter 
partial  sequences  for  LPI-2  that  gave  better  than  AWGN  performance  in  the  Harbor 
scenario.  Thus,  the  AWGN  analysis  and  the  ELB  analysis  provide  benchmarks  for 
assessing  the  LPI  multipath  performance.  This  stopping  point  is  the  natural  starting 
point  for  Phase  4.  Specifically,  Phase  4  should  undertake  the  following:  (1)  extend 
the  analysis  to  additional  LPI  systems,  (2)  search  for  the  presence  of  a  noise  floor,  (3) 
increase  the  bit  rate,  (4)  extend  the  ELB  analysis  to  encompass  the  RF  environmental 
effects,  (5)  assess  the  LPI  performance  in  multipath,  (6)  develop  these  results  into 
credible  RF  mission-planning  tools. 

1.2  Notation 

The  notation  was  selected  to  be  consistent  with  as  much  of  the  literature  as  possible 
and  to  permit  the  reader  to  “drop  in”  on  an  equation  without  a  lot  of  notational 
baggage.  We  follow  Staley’s  notation  in  his  excellent  thesis  [38]  wherein  the  subscripts 
“2”’  and  “i?”  denote  the  transmitter  and  receiver.  Table  1.1  was  designed  with  these 
rules  and  with  an  effort  to  reduce  notational  collisions. 
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Table  1.1:  Notation. 


Symbol 

Description 

/ 

frequency  (Hertz) 

fc 

carrier  frequency  (Hertz) 

u 

fade  rate  (Hertz) 

fn 

maximzJ  Doppler  shift  (Hertz) 

A/d 

relative  Doppler  shift  (Hertz) 

G'P 

gain  of  transmitter’s  antenna 

Gr 

gain  of  receiver’s  antenna 

Go 

free-spaoe  propagation  gain 

hx 

transmitter’s  antenna  height  (feet) 

hR 

receiver’s  antenna  height  (feet) 

h(t,T) 

Input  Delay-Sprezwi  function 

hRF(t,r) 

RF  Input  Delay-Spread  function 

hs  {m,k) 

Baseband  Digital  Input  Delay-Spread  function 

Os 

oversampling  factor 

A 

wave  length  (meters) 

PfU) 

Doppler  scattering  function 

Pt{t) 

power  delay  profile 

Pxit) 

transmitter’s  pulse-shaping  filter 

Pr(^) 

receiver’s  matching  filter 

Pr 

power  transmitted  (Watts) 

Pr 

power  received  (Watts) 

Pq 

power  in  fi:ee-space  (Watts) 

■Pxx(/) 

power  spectrum  of  {ir(t)} 

Rxxitiyt2) 

covariance  of  {x(t)} 

R 

information  bit  rate  (bits/sec) 

t 

time  (seconds) 

T 

symbol  duration  (seconds) 

Tc 

chip  duration  (seconds) 

r 

delay  (seconds) 

At 

relative  delay  (seconds) 

e. 

grazing  angle  (radians) 

'S'xx(/1j  /2) 

two-dimensional  spectrum  of  {x(t)} 

mr) 

Doppler-Delay  function 

W 

S-d6  bandwidth  (Hertz) 

®RCW 

raised  cosine  filter 

The  Fourier  transform  of  a  signal  s{t)  is 

s(/)  =  r 

oo 

The  convolution  of  s{t)  with  a  filter  h{t)  is 

h  *  s(t)  =  f  h{t  —  T)s{r)dT. 

J—OO 

A  random  process  is  denoted  by  {a:(t)}.  The  covariance  of  {x(t)}  is 

—  .^[3:(ti)3:(t2)]) 

where  the  complex  conjugate  is  denoted  by  the  overline.  If  {x(t)}  is  harmonizable, 
then  it  admits  a  representation  by  a  stochastic  Fourier  measure: 

/OO 

e+j2^ftdX{f). 

'OO 
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Equivalently,  i?xa:(ii>  ^2)  admits  a  representation  as  the  two-dimensional  Fourier  trans¬ 
form  of  the  covariance  of  X{f): 

R^ih,  h)  =  r  r  ,  h), 

J— CO  J~oo 

where  Sxxifiifi)  —  ^[-^(/i)'^(/2)]-  If  is  wide-sense  stationary  (WSS),  then 

it  is  standard  to  write  Rxxituh)  ~  Rxx{ti  —  h)-  This  is  equivalent  that  the  two- 
dimensional  spectrum  is  supported  on  the  main  diagonal  in  the  frequency  plane: 
Sxx  =  Pxxifi)  ^(/i  -  /z))  where  Pxx{f)  is  the  usual  power  spectrum. 

Matrix  and  vector  notation  is  standard.  A  vector  x  is  in  column  form  with  its  size 
denoted  x  ~  A*  x  1,  where  Nx  denotes  the  number  of  rows  of  x.  Matrices  are  upper 
case  A  with  dimensions  indicated  as  A  ~  M  x  iV.  The  matrix  transpose  is  and  the 
conjugate  transpose  is  A^.  The  ubiquitous  MATLAB  grammar  is  also  used.  A(:,n) 
denotes  the  nth  column  of  matrix  A.  A{m, :)  denotes  the  mth  row.  The  index  vector 
is  denoted  [1  :  A]  =  [12  ...  N].  Then  x(3  +  [1  :  N])  =  [a:(3)  r(4) . . .  x{N  4-  3)]^. 
The  “floor”  function  rounds  a  real  number  toward  —00:  The  “argmax”  function  is 
the  inverse  of  the  max  function.  If  h  :  R  and  h^ax  =  max{h(u;)  :  a;  G  0}  then 

argmax[h]  :=  {w  G  :  h(w)  =  hmax}- 


1.3  The  Baseband  System 

Johnson’s  actual  LPI  system  employs  several  non-trivial  techniques  to  hide  the  fea¬ 
tures  of  the  radiated  power  [19].  His  system  cleverly  compensates  for  these  parameter 
variations.  By  assuming  optimal  receiver  operations,  we  ignore  these  variations  and 
work  from  a  fixed  configuration.  Thus,  we  estimate  LPI  performance  for  the  generic 
baseband  transmitter  and  receiver  shown  in  Figure  1.1.  The  baseband  assumptions 
(B-1,  B-2,  . . .  B-5)  associated  with  this  figure  are  listed  next: 


ti(») 


Figme  1.1:  Baseband  transmitter  and  receiver  model 
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B-1  Linear  Modulation.  The  signal  SRp(t)  transmitted  at  RF  is  obtained  from  the 
baseband  signal  s{t)  as 

SRF{t)  = 

B-2  Linear  Demodulation.  The  received  baseband  signal  e{t)  is  demodulated  from 
the  RF  signal  as 


e(t)  =  e  /  /iRp(t,r)sRF(t  -  r)(ir, 

J—OO 

where  /iRF(t,r)  denotes  the  RF  Input  Delay-Spread  function. 

B-3  Perfect  Synchronization.  The  receiver  and  transmitter  operate  from  a  master 
clock. 

B-4  Matched  Filter.  The  convolution  of  the  transmitter  filter  pr(0  and  the  re¬ 
ceiver  filter  pR(t)  is  the  raised  cosine  filter  [34,  pages  63-65]:  PR*PT{t)  = 

B-5  Slow  Fading.  The  Baseband  Input  Delay-Spread  function  has  the  form: 

N 

n=l 

where  the  {a„(t)}’s  are  JWSS  complex  Gaussian  with  mean  and  spectral  shape 
determined  by  the  propagation  environment  and  the  RF  modulations.  The 
“slow-fading”  assumption  is  that  the  maximal  spectral  extent  of  the  {an(t)}’s 
is  “small”  with  respect  to  the  bandwidth  of  the  receiver  filter  pR{t). 

Following  [30]  and  [18],  the  baseband  channel  input  has  the  complex  baseband  form: 

OO 

sW=  s  ^kPrit-kTc), 

k=—oo 

where  the  s^’s  will  be  the  “chips”  (sj,  =  ±1)  from  various  PN  sequences.  The  base¬ 
band  receiver  output  is  modeled  as 

yii)  =Pr*  e{t)  +Pr*  ri{t), 

where  is  a  complex-valued,  zero-mean,  white  Gaussian  random  process  with 

two-sided  power  spectral  density  2No.  The  slow  fading  assumption  permits  the  ap¬ 
proximation  [1,  Eq.  3.2]: 

OO  N 

yit)^  E  Ylskanit)xRcit-Tr,-kT,)  -f-  PR*ri{t).  (1.1) 

k=—oo  n=l 
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For  the  additive  white-noise  Gaussian  channel,  r)  =  a5{t).  In  this  case,  the 
receiver’s  output  is 


y{t)  «  a  ^  -  kTc)  +  Pr*  v{t). 

k=—oo 

By  B-3,  ideal  symbol  timing  is  possible.  That  is,  the  receiver  may  sample  y{t)  at 
instants  tm  =  rnTjOsj  where  Os  denotes  the  oversampling  factor.  With  0*  =  1, 
digital  input  to  the  demodulator  has  no  ISI  and  then  has  the  form: 

y(^m)  ~  ”1“  9m'>  (^•^) 

where  {pm}  is  Gaussian  IID. 
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LPI  Performance  in  AWGN 


This  section  establishes  a  basehne  performance  of  the  LPI  systems.  A  generic  base¬ 
band  transmitter  and  receiver  are  assumed  (Section  1.3).  The  channel  is  additive 
white  Gaussian  noise  (AWGN).  For  each  noise  level  and  each  receiver  parameter,  we 
obtain  an  LPI  performance  vector  [Pe  E^/Nq  R/W  SNR]^.  The  SNR  is  used  as  a  sim¬ 
ple  measure  of  covertness.  As  we  sweep  over  the  noise  levels  and  detector  parameters, 
the  resulting  performance  vectors  generate  an  LPI  performance  manifold.  Slices  of 
this  manifold  provide  handy  graphics  of  the  LPI  performance.  For  example,  we  wiU 
see  that  each  LPI  system  determines  a  two-dimensional  manifold  in  Pe  x  R/W  x  SNR 
space.  If  a  system  engineer  or  RF  mission  planner  needs  an  operational  SNR  and 
probability  of  symbol  error,  then  the  spectral  efficiency  is  not  also  independent  but 
located  on  the  manifold  at  a  point  determined  by  the  SNR  and  Pg.  With  the  spectral 
efficiency  then  fixed,  a  bit-rate  requirement  will  set  the  actual  bandwidth.  That  is, 
the  LPI  geometry  determines  RF  and  receiver  pmameters.  Thus,  a  useful  LPI  mission 
tool  is  developed  for  each  of  the  LPI  receivers.  A  glance  at  the  table  of  contents  will 
show  that  for  each  LPI  system,  we  estimate  the  theoretical  performance,  examine  the 
LPI  settings  in  noise,  and  collect  all  these  results  in  the  final  system  performance  (the 
LPI  manifold). 

2.1  Parameterizing  the  LPI  Systems:  [Np^N^^Ny\ 

The  generic  LPI  scheme  is  a  variant  of  non-coherent  M-ary  orthogonal  signaling  using 
symbol  vectors  ui,  ua,  .. . ,  ui.  The  u/s  are  obtained  from  a  generalized  partition 
of  a  PN  sequence  u  consisting  of  ±l’s  of  length  Currently,  the  u^’s  aU  have  the 
same  length  Np.  Np  is  also  called  the  partial  length  and  1  Np  <  N^.  The  indices 
of  adjacent  symbol  vectors  may  overlap  by  length  0  <  Na  <  Np.  For  1  =  1,  L 

Ui  =  u(  (1  -  1)  X  (Np  -  Na)  -b  [1  :  Np] ), 
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where  cyclic  wrapping  of  u  is  assumed  when  necessary.  Although  the  symbol  vectors 
are  not  orthogonal,  the  goal  is  to  have  small  cross  correlations.  This  tends  to  be  true 
for  large  Np  and  because  of  the  “white  noise”  nature  of  the  PN  sequences.  Thus, 
the  variables  {Np,N£^,Ny,)  can  be  viewed  as  three  “knobs”  that  dial  in  the  desired 
performance  from  the  corresponding  LPI  receiver. 

This  report  typically  fixes  at  =  2^®  but  varies  the  partial  length  Np  and  overlap 
Na-  LPI-1  takes  the  overlap  to  be  Na  =  0  so  that  the  symbol  vectors  are  a  partition 
of  the  PN  sequence  u.  In  contrast,  LPI-2  uses  maximal  overlap  Na  —  Np  —  1  so 
that  many  more  symbol  vectors  are  available  in  comparison  to  LPI-1.  We  will  see 
that  this  gives  LPI-2  a  higher  spectral  efficiency  but  at  the  expense  of  covertness  and 
sensitivity  to  multipath. 


2.2  LPI-1:  An  M-ary  LPI  Scheme  with  =  0 

LPI-1  is  a  variant  of  non-coherent  M-ary  orthogonal  signaling  using  symbol  vectors 
ui,  U2,  . . . ,  U£,.  The  uj’s  are  obtained  from  a  partition  of  a  PN  sequence  u.  Conse¬ 
quently,  the  symbol  vectors  are  not  orthogonal  but  have  small  cross  correlation.  The 
baseband  assumptions,  with  B-3  modified  to  include  the  knowledge  of  these  symbol 
vectors  and  their  transmission  times,  give  the  vector  form  of  the  detector’s  input  (See 
Equation  1.2): 

y  =  [y(im)]  =  auz  +  g. 

For  the  simulations,  we  take  u^’s  to  be  non-overlapping  with  length  Np.  If  the  entire 
PN  sequence  has  length  N^  then  there  are 

L  =  Qoov{Nu/Np) 

symbol  vectors.  This  L  will  be  used  throughout  the  LPI-1  analysis.  If  NpL  <  N^ 
then  the  trailing  chips  are  dropped.  Form  the  matrix 

17  =  [Ui  U2  ...  Ui]  rsj  NpXL. 


The  detector  output  may  be  written  as 

T=  argmax(|17-*^y|). 


2.2.1  Pe  Estimates 

Because  the  detector  uses  magnitude  only,  the  channel  loss  may  be  scaled  into  the 
noise.  Let  y  denote  this  scaled  input  vector  to  the  detector.  Then 


y  =  Ui  +  g, 
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where  g  ~  Af{0,agl).  The  SNR  at  the  detector’s  input  is  the  ratio  of  mean  energy 
in  the  s5Tnbol  vector  to  the  noise  energy: 

SNR  =  4- 

We  estimate  the  probabihty  of  error  as  a  function  of  this  SNR.  Let  r  denote  a  scaled 
version  of  the  detected  vector: 


Nr,  ^  Nr,  ^  Nr,  ^ 


J-^p  J.'lp 

Following  Miller  [27,  page  77],  the  probabihty  density  of  r  given  I  is 

1 


=  7r«,det|fl|^P  -  r,))  . 

1 


where 


Ti  =  —Uui  and  R 


N2 


Assuming  all  L  symbols  equally  hkely  gives  the  probabihty  of  symbol  error  as 


^  1=1 


1=1 


The  problem  is  to  estimate  P[l  ^  Z|Z]  to  determine  the  symbol  error  Pg. 

An  Orthogonal  Approximation-.  Because  U  NpX  Lis  formed  from  a  PN  sequence, 
it  is  approximately  orthogonal:  U^U  ^  NpXl^.  Then,  assuming  ti  was  transmitted, 
r  admits  the  approximate  distribution: 


r  « 


1 


l+AfiO^g/Np) 


l'  =  l 

- 


That  is,  the  components  of  r  are  approximately  independent  and  complex  Gaussian. 
Then  the  components  of  |r|  are  approximately  independent  and  Ricean  [34,  page  47]: 


Pv\i{r)  «  l[0,co)(r) 


r 


g-(P+l)/2cr2  P  ^  I 

e-rV2.^  r  ^  I 


where  cr^  =  crg/{2Np).  Then  the  joint  probabihty  density  of  |r|  is 


Pl,2,...,L|i(»'l>  ^2,...,  tl)  «  Pi\iiri)p2\i{r2) . .  .PL\l{rL)- 

Then  the  probabihty  of  a  correct  decision  P\l  =  /[/]  is  the  integral  of  the  joint  density 
over  the  set 

72.(Z)  =  {r  6  :  rj  >  rii  for  /'  =  1, . . .  X,  Z'  ^  Z}. 
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More  formally  [37,  page  179], 


P[T  =  l\l]  =  Prob[|r(/)|  >  |r(0|  for  =  1, . . . ,  L,  /'  ^  /] 
=  /  .  ^2,  •  •  • ,  ^l)  <iridr2 . . .  dr^ 

77^(0 

^  rr 

J  dripiii{ri)  JJ  J  drv  pi.ii{r[) 

vil 

=  jT  /o(r/cr^){l  - 


=  ^  V  / 


!72<r=(I'+l) 


This  result  shows  that  the  probability  of  a  correct  decision  is,  with  respect  to  this 
approximation,  not  dependent  on  1.  Then,  for  any  I  there  holds 


=  7  'I'l  =  -PpT^  iW  =  1  -  /’P=  i|i|. 

l'=l 


The  following  simulations  demonstrate  the  quality  of  this  approximation  using  symbol 
vectors  on  the  order  of  Johnson’s  systems. 
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2.2.2  Pe  for  [Np,  Na,  N^]  =  [1024, 0, 2048] 

A  low-dimensional  simulation  illustrates  this  PPM  system.  The  PN  sequence^  u  has 
length  iV„=2048.  With  partial  length  Np  =  1024,  the  correlator  is  nearly  orthogonal: 

1.0000  -0.0215  ■ 

-0.0215  1.0000  ■ 

Figure  2.1  compares  the  simulation  with  the  Pg  estimate.  An  excellent  fit  is  obtained 
provided  we  observe  that  the  number  of  simulations  (100,000)  limit  the  trustworthy 
estimates  to  Pg  >  10"“*.  A  rule-of-thumb  is  that  U  becomes  “orthogonal”  provided 
Np  >  iVu/10  [1].  For  system  applications,  increasing  Np  permits  a  lower  SNR  but 
reduces  the  bit  rate.  These  trade-offs  are  explored  next. 


PN(1 1,5001);  PartiBl  Sequence  Length=1024;  L=2 


rTTT^ - 1 - r - 1 - j - 1  I - 1 - j 

O  Simulation 

o  1 

- Approximation  ' 

“O 

r  *o 

I 

o 

r  o 

- 

b. 

- « - 1 - 1 - > - 1 - 1 _ 1 _ 

0\  r 

1 


-25  -24  -23 


-22  -21  -20  -19  -18 

SNR;  Number  of  simulations^ 100000 


-17  -16  -15 


Figure  2.1:  Comparison  of  Pg  (symbol  error)  simulations  and  approximation  for 
[Np,Na,Nu]  =  [1024,0,2048]. 


^Generated  as  u  =  pn_scr(ll,  5001,  1)  from  [31] 
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2.2.3  P.  for  [Afp,  Na,  Af„]  =  [1024, 0, 2i«] 

A  simulation  more  in  line  with  one  of  Johnson’s  LPI  systems  uses  a  PN  sequence^  u 
of  length  Nu  =  2^®.  With  partial  length  Np  =  1024,  we  are  well  below  the  10%  rule. 
Consequently,  the  approximation  is  more  crude  with 

j^U-U  =  I^E. 

where  most  of  the  elements  of  E  are  still  small  (mean  magnitude  «  0.02).  However, 
a  few  elements  have  magnitude  fa  0.1.  Although  the  orthogonal  approximation  is 
degraded,  it  still  gives  good  Pg  estimates  as  shown  in  Figure  2.2.  Although  LPI-1 
with  L  =  64  has  6  times  the  bit  rate  as  the  system  with  L  =  2,  the  SNR  must  increase 
for  comparable  symbol  error  rates,  as  can  be  seen  by  comparing  Figures  2.1  and  2.2. 


Figure  2.2:  Comparison  of  Pg  (symbol  error)  simulations  and  the  approximation  for 
[iVp,iVA,iV„]  =  [1024,0,21®]. 


^Generated  as  u  =  pnjscr(16, 210013,  iVu,!)  from  [31]. 
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2.2.4  System  Performance:  [Np,NA,Nu]  =  [iVp,0, 2^®] 

This  section  presents  performance  trade-offs  for  the  LPI-1  system.  We  start  by  com¬ 
paring  several  partial  lengths  Np  in  Figure  2.3.  The  longer  partial  lengths  permit 
operation  in  lower  SNR  but  reduce  the  information  bit  rate  R. 


Figure  2.3:  Comparison  of  Pg  (approximations)  against  SNR  for  partial  lengths  Np 
for  [N^,N^]  =  [0,2i«]. 


To  quantify  this  statement,  let  each  chip  have  duration  Tc.  With  the  raised-cosine 
pulse,  the  3-dB  bandwidth  is  W  =  1/Tc.  Each  symbol  vector  has  duration  T  —  NpTc 
so  that  the  information  bit  rate  is 


_  hg^jL) 

NpTa 


[bits/sec . 


The  spectral  efficiency  is  [34,  pages  282-284] 


Figure  2.4  plots  the  resulting  LPI  performance  curves  as  a  function  of  the  partial 
lengths  Np  This  plot  is  really  a  “skeleton”  of  the  LPI  performance  surface.  That  is, 
LPI  performance  variables  SNR,  Pe,  R/W  are  constrained  to  lie  a  surface.  For  clarity, 
we  only  plot  the  associated  Hnes.  This  surface  reveals  the  LPI  the  trade-offs:  more 
covertness  or  low  SNR  requires  longer  partial  lengths  that  reduce  R.  This  surface 
is  also  assists  RF  mission  planning.  For  example,  a  glance  at  Figure  2.4  shows  that 
R/W  is  a  function  of  SNR  and  Pe.  If  the  covertness  requirements  set  the  SNR  near 
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-25  dB  and  the  error  requirements  are  10~^  then  R/W  «  0.1  bits/sec/Hz.  To  get 
2400  bits  per  second  dictates  a  bandwidth  of  approximately  24  kHz.  This  surface  also 
shows  the  RF  mission  planner  what  trade-offs  can  be  made  if  the  desired  operational 
point  is  not  feasible.  For  example,  if  R/W  must  be  near  1  bits/sec/Hz  then  we  can 
either  (1)  relax  to  10”^,  (2)  relax  SNR  to  -20  dB,  (3)  or  simultaneously  relax 
to  10"*  and  the  SNR  to  -23  dB. 

Estimaled  LPI  Pafonoance:  PN  )engt}>=^S36 


Figure  2.4:  LPI  performance  as  function  of  partial  lengths  Np  for  [N^,  iV,,]  =  [0,2^®]. 


We  close  with  the  standard  Ei/No  plots.  E\,/Nq  is  derived  from  the  SNR  as 


Eb  ^  Np 
No  ~  log2(L) 


X  SNR. 


Figure  2.5  plots  Pe  as  a  fimction  of  Ei/Nq.  This  figure  agrees  with  Figure  3.34  in 
Sklar  [37]  for  non-coherently  detected  M-ary  orthogonal  signaling,  noting  the  bit  error 
versus  symbol  error  difference  between  the  plots.  This  performance  curve  apphes  to  all 
M-ary  orthogonal  signaling,  regardless  of  spectral  efficiency.  Figmre  2.6  plots  spectral 
efficiency  against  Eh/No.  Comparing  Figure  2.6  with  the  bandwidth-efficiency  plane 
in  Figure  7.6  of  [37],  we  see  that  LPI-1  performance  exists  in  the  swath  of  Eb/No 
firom  7  to  11  dB,  which  agrees  with  MFSK  performance.  However,  the  LPI  nature 
of  the  waveform  comes  at  the  expense  of  spectral  efficiency,  which  is  much  less  than 
that  of  the  corresponding  MFSK  waveforms.  Also,  all  of  the  LPI-1  waveforms  in 
Figure  2.5  and  2.6  have  the  same  bandwidth,  in  contrast  to  MFSK  in  Figure  7.6  of 
[37].  So,  as  L  decreases  in  LPI-1,  the  spectral  efficiency  also  decreases.  These  figures 
reveal  a  trade-off.  Although  the  longer  sequences  are  more  covert  in  terms  of  the 
received  SNR,  the  lower  data  rates  associated  with  these  sequences  result  in  a  loss  in 
performance  with  respect  to  Eb/No  and  spectral  efficiency. 
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Figure  2.5:  Comparison  of  Pj.  (bit  error)  against  Ei/Nq  for  partial  lengths  Np.  — 
2^®.  Each  line  from  left  to  right  is  determined  by  the  partial  lengths  7^^=1024,  2048, 
4096,  8192,  16384. 


10‘ 


Euimated  LPI  Berfotnunce  at  P 


OiOM 

03M< 

i  1 - 

O40» 

0<in 

. - . -.J _ L. 

9.5  9  9.5 

Ej/Nq  (dB);  PN  lcn|th=«5536 


Figure  2.6:  Comparison  of  R/W  against  Eb/No  for  partial  lengths  Np  for  TVa  =  0 
and  iV„  =  2i«. 
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2.3  LPI-2:  A  CCSK  LPI  Scheme  with  N/^  =  Np  —  1 

Johnson  invented  an  LPI  system  based  on  cyclic  code-shift  keying  (CCSK)  [19].  As 
before,  we  assume  Johnson’s  receiver  compensates  for  his  LPI  parameter  variations. 
The  CCSK  scheme  starts  from  a  given  PN  sequence  u  of  length  JV^.  The  sjunbols 
1  =  0,  1,  .,.Nu  — I  are  mapped  to  cyclic  shifts  of  u.  Thus,  the  Zth  symbol  vector  for 
CCSK  is 

uj  =  tt(rem(l  +  [1 :  A^],  Au)) 

with  L  =  Nu  symbol  vectors.^  To  increase  the  bit  rate,  Johnson  invented  a  partial 
CCSK  system  by  exploiting  the  almost  orthogonality  of  u  by  only  transmitting  the 
first  Ap  chips  of  u^.  That  is,  the  Zth  symbol  vector  for  partial  CCSK  is 

ui  =  u(rem(i  •+•  [1 ;  Ap],  A^)). 

These  L  symbol  vectors  each  carry  log2(A„)  bits  of  information.  The  baseband  as¬ 
sumptions  give  that  the  vector  form  of  the  detector’s  input  (See  Equation  1.2): 

y  =  [y{tm)]  =  aui  +  g. 


The  detector  output  is 

T  =  argmax(|u  *c  y|)) 
where  the  subscript  “c”  denotes  circular  correlation. 

2.3.1  Pe  Estimates 

Because  only  the  magnitude  matters  for  the  detector,  we  scale  the  channel  loss  into 
the  noise.  If  y  denotes  this  scaled  input  vector  to  the  detector  then 

y  =  uz  +  g, 

where  g  ~  A^(0,  CTgl).  The  SNR  at  the  detector’s  input  is  the  ratio  of  the  mean  energy 
in  the  symbol  vector  to  the  variance  of  the  complex  additive  Gaussian  noise: 

SNR  =  4- 

We  estimate  the  probability  of  error  as  a  function  of  this  SNR.  Introduce  the  A«  x  Np 
circulant  matrix  Cu  to  model  the  cychc  correlation  with  u: 

U  *c  Uz  =  CuUj. 

®  “rem”  denotes  the  remainder  function.  For  brevity,  indexing  starts  at  0. 
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Let  r  denote  a  scaled  version  of  the  detected  vector: 


The  scaling  is  rigged  so  that  the  maximum  of  the  correlation  is  1.  Following  Miller 
[27,  page  77],  the  probability  density  of  r  given  I  is 


p(r|/) 


•exp  (-(r  -  r,)^i7  ^(r  -  rj)) 


TT-^p  det[i7] 

where  rj  =  C^ui/Np  and  the  noise  has  the  covariance  matrix 

Then  r  is  approximately  distributed  as 

f  l+^-(0,,7“/JVj)  l'  =  l 

"  ~  \  Af{0,Tl/N^)  I'jil- 

That  is,  the  components  of  r  are  approximately  independent  and  complex  Gaussian. 
This  is  the  approximation  used  in  Section  2.2.1.  With  each  symbol  equally  likely, 

I  Nu-l 

p>  =  w^ 

t'=0 

The  probability  of  a  correct  decision  is  approximately 

P[T=l\l]  « 


iV^-1 


=  E 


(-1)^'  /  iVu  - 1 

i'  +  i[  v 


-r/2cr2(/'+l) 


where  =  cr^/(2iV^).  This  shows  the  probability  of  a  correct  decision  is,  with  respect 
to  this  approximation,  not  dependent  on  1.  Then,  for  any  I, 

p^=p[r5^/ii]  =  i_p[r=i|/]. 

The  quality  of  this  approximation  is  assessed  in  the  following  subsections. 


2.3.2  Pe  for  [iVp,  iVA,  iVj  =  [1024, 1023, 2^®] 

The  PN  sequence^  u  has  length  iVu=65,536.  The  partial  lengths  are  set  to  iVp— 1024 
for  comparison  with  the  LPI-1  (M-ary)  simulations.  In  the  LPI-1  simulation  of  Sec¬ 
tion  2.2.3,  each  iVp=1024  chips  carried  6  bits.  In  this  LPT2  simulation,  each  partial 

‘‘Generated  by  u  =  pn_scr(16, 210013,  iVu,  1)  from  [31]. 
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sequence  of  iVp=1024  chips  carries  16  bits  of  information.  Consequently,  the  bit  rate 
for  LPI-2  has  increased  by  a  factor  of  2.7=16/6  over  LPI-1.  In  comparison  to  standard 
CCSK,  the  bit  rate  for  LPI-2  has  increased  by  a  factor  of  64.  However,  increasing 
the  bit  rate  incurs  a  loss  of  covertness.  Figure  2.7  illustrates  this  point  by  comparing 
LPI-2  against  LPI-1.  Although  the  partial  sequence  length  Np  is  the  same  for  both 
LPI-1  and  LPI-2,  there  are  2^°  more  s3anbol  vectors  in  the  LPI-2  symbol  set.  There¬ 
fore,  more  potential  symbol  errors  are  possible  in  LPI-2  versus  LPI-1.  As  a  result, 
LPI-2  requires  a  higher  SNR  for  comparable  Pg  performance. 

LPl-l=[l024,0.6553ej;U»I-2=[l024.l023.6SS36J 


Figure  2.7:  Comparing  Pg  for  LPI-1  (M-ary)  [Ap,  JVa,  iVu]  =  [1024,0,2^®]  against 
LPI-2  (partial  CCSK)  with  [Ap,  Aa,  A„]  =  [1024, 1023, 2^®]. 
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2.3.3  Pi  for  [jVp,  Na,  N^]  =  [2i«,  2“  -  1, 2'S] 

Although  a  special  case  of  partial  CCSK  with  Np  =  N^,  Figure  2.8  compares  the 
theoretical  approximations  with  measurements  obtained  by  George  B.  Johnson  in 
November  1999.  Agreement  between  theory  and  measurements  demonstrate  that 
reasonable  modeling  of  the  COTS  LPI  system  has  been  obtained.  Here  the  bit  error 
Pfr  is  plotted  and  relates  to  the  symbol  error  Pe  as  [34,  page  311]: 

2^-1 

where  M  =  log2(L)  is  the  number  of  bits  carried  by  each  symbol  vector.  What  is 
striking  about  this  figure  is  the  low  SNR  that  may  be  obtained  using  the  long  PN 
sequences.  This  reduction  in  SNR  costs  in  terms  of  the  bit  rate.  Consequently,  the 
next  section  completes  this  examination  of  CCSK  by  making  the  bit  rate  versus  SNR 
trade-off. 


Figure  2.8:  Compeiring  bit  error  estimates  between  the  measurements  and  theory  for 
standard  CCSK:  [Np,  Na,K]  =  [2^®, 2^®  -  1,2^% 
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2.3.4  System  Performance:  [Np,Na,Nu]  =  [Np,Np  —  1,2^®] 

This  section  offers  performance  trade-offs  for  the  LPI-2  (partial  CCSK)  system.  For  a 
fixed  PN  length  iVit=65536,  we  compare  several  partial  lengths  Np  in  Figure  2.9.  The 
longer  partial  lengths  permit  operation  in  lower  SNR,  but  reduce  the  information 
bit  rate  R.  To  quantify  this  statement,  let  each  chip  have  duration  Tc.  With  the 
raised-cosine  pulse,  the  3-dB  bandwidth  is  W  =  l/Tc-  Then  each  symbol  vector  has 
dmration  T  =  NpTc.  Recall  that  each  symbol  vector  carries  log2(iVu)=16  bits.  The 
information  bit  rate  is 


R 


[bits/sec]. 


The  spectral  efficiency  is  [34,  pages  282-284]: 


W 


log2(-^u) 

Nr, 


[bits/sec/Hzj. 


Figure  2.10  plots  the  resulting  LPI  performance  curves  as  a  function  of  the  partial 
lengths  Np.  The  more  covert  operation  (SNR'd)  requires  longer  partial  lengths  with 
the  subsequent  reduction  of  R.  Eb/No  is  given  by 


No 


Np 

l0g2(iV„) 


X  SNR. 


Figure  2.11  plots  Pe  as  a  function  of  Eb/No.  This  plot  reveals  how  misleading  Eb/No 
performance  measures  are  for  an  LPI  system.  Although  the  longer  sequences  operate 
in  lower  SNR  environments  as  Figme  2.9  shows,  there  is  a  decrease  in  the  bit  rate 
such  that  Eb/No  is  identical  for  all  these  sequences.  Figure  2.11  also  does  not  have  the 
familiar  form  of  Figure  3.34  in  [37]  for  non-coherently  detected  orthogonal  signaling 
because  each  of  the  partial-length  symbols  LPI-2  encodes  the  same  number  of  bits. 
Each  curve  in  Figure  2.11  falls  on  the  M  =  2^®  curve  in  Figure  3.34  in  [37],  although 
the  spectral  efficiently  differs  for  each  waveform  as  Figure  2.10  shows. 
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Figure  2.9:  Comparison  Pg  (symbol  error)  against  SNR  for  partial  lengths  Np: 
[Np,NA,N^]  =  [Np,Np-l,2^% 


Estimaied  LPl  ^rfonnancc  PN  lenglh:=6SS36 


Figure  2.10:  LPI-2  performance  as  function  of  partial  lengths  Np-.  [Np,N^,Nu]  = 
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Eslimatul  LPI  Perfomiance:  PS  length=6SS36 


Figure  2.11:  Comparison  of  Pg  (sjunbol  error)  against  Eb/No  for  partial  lengths  Np: 

[Np,NA,N^]  =  [Np,Np-l,2^% 


2.4  Parameterizing  the  LPI  Performance 

This  section  concludes  by  showing  the  conjectured  mapping  of  the  LPI  parameters 
[iVp,  TV^jiVu]  to  the  system  performance  [Pe,i?/W,  SNR]  in  additive  white  Gaussian 
noise  with  variance  (Xg.  The  approximations  are  good  for  LPI-1  and  LPI-2  over 
the  [ATpjiVAj-Nu]  values  used  in  this  section.  It  remains  to  be  seen  how  well  the 
approximations  hold  over  the  entire  three-parameter  LPI  family.  Given  [Np,  JVa,  Nu] 
assume  that  L  symbol  vectors  ui,  ...,  U£  are  transmitted.  The  computation  of 
L  depends  on  the  generalized  partition  and  the  cyclic  wrapping  of  u.  If  no  cyclic 
wrapping  is  used  then  L  is  bounded  as 

-  Np-N^ 


For  LPI-1,  there  are  L  =  floor (iVu/iNTp)  symbol  vectors.  For  the  CCSK  in  LPI-2,  cyclic 
wrapping  gives  L  =  Ny,  symbol  vectors. 

The  SNR  is  the  ratio  of  the  mean  energy  in  the  symbol  vector  to  the  variance  of 
the  complex  additive  Gaussian  noise: 


SNR  = 
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The  spectral  efficiency  is  the  rato  of  the  bit  rate  R  to  the  bandwidth  W: 


^  ^  logzCi^) 

W  Np  ’ 


where  each  chip  has  duration  Tc-  For  the  probability  of  symbol  error,  the  current 
conjecture  is  that 


Pe 


L-1  ( 

i-Y'  - - 


1)^' 


This  approximation  appears  to  be  good  for  LPI-1  with  L  =  floor(iVu/JVp)  and  LPI-2 
with  L  —  Nu.  Assessing  the  range  of  validity  of  these  formulas  is  one  the  Phase  4 
tasks  mentioned  in  the  beginning  of  this  report.  Likewise,  extending  these  trade-off 
formulas  to  multipath  environments  is  another  Phase  4  task.  Finally,  a  synthetic  case 
of  multipath  is  to  overlap  the  symbol  vectors  in  time.  This  increases  the  bit  rate  but 
should  also  increase  the  SNR  and  error  rate.  Quantifying  the  trade-off  between  bit 
rate  and  covertness  is  a  Phase  4  task. 
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The  ELB  Channel  Models 


Section  2  determined  LPI  performance  in  additive  white  Gaussian  noise.  For  real- 
world  operations,  the  LPI  performance  must  be  assessed  in  operationally  credible 
multipath.  Referring  to  the  baseband  system  in  Figure  1.1,  this  means  that  an  Input 
Delay-Spread  function  h{t,  r)  needs  to  be  specified,  in  addition  to  the  additive  noise. 
Extended  Littoral  Battlefield  (ELB)  multipath  provides  an  excellent  “testbed”  for 
such  an  assessment.  This  section  presents  three  ELB  channel  models  to  be  used  in 
the  multipath  simulations  of  Section  4: 

Bluewater.  Mobile  ship  in  bluewater  to  shore  site. 

Hsirbor.  Mobile  ship  in  urban  harbor  to  shore  site 

HMMWV.  Mobile  HMMWV  in  urban  environment  to  shore  site. 

These  ELB  channel  models  are  derived  from  the  Naval  Communication  Channels 
(NCC)  Dataset.  The  NCC  Dataset  contains  over-the-air  measurements  of  a  mobile 
HMMWV  and  a  mobile  Navy  ship  transmitting  to  a  fixed  shore  site  at  1  Mbps  in 
the  225  to  400  MHz  military  UHF  band  collected  at  San  Diego,  Cahfornia  during 
January  and  February  of  1998  [1].  As  such,  these  experiments  are  representative  of 
urban  Marine  Corps  and  the  Naval  communications  in  the  ELB. 

This  section  analyzes  this  ELB  dataset.  The  immediate  goal  provides  collection 
of  credible  ELB  channel  models  for  LPI  multipath.  The  long-term  goal  develops  an 
understanding  of  ELB  propagation  for  RF  mission  planning  and  LPI  deployments. 
Section  3.1  reviews  general  channel  models.  Section  3.2  specializes  these  models  to 
the  ELB  propagation  environment  Section  3.3  describes  the  receiver  and  transmitter, 
the  data  acquisition  process,  the  routes  of  the  ship  and  HMMWV.  Comparison  with 
measurements  indicate  that  a  Ricean  model  provides  an  upper  bound  on  the  received 
power  levels.  Section  3.4  describes  channel  estimation  from  noisy,  distorted,  digitized 
versions  of  the  RF  signal  that  probed  the  ELB  channels.  The  remaining  sections  each 
describe  the  bluewater,  harbor,  and  HMMWV  channel  models.  The  goal  is  to  reduce 
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each  channel  to  a  table  listing  each  path,  its  power  level,  delay,  Doppler,  fade  rate, 
and  fading  processes.  These  tables  are  the  end  product  of  this  effort.  The  channel 
modeler  can  simply  turn  to  these  tables  to  program  a  channel  emulator  or  continue 
to  the  next  section  that  puts  the  LPI  systems  into  these  channels. 

3.1  General  Channel  Models 

Bello’s  seminal  1963  paper  estabhshed  the  nomenclature  and  analysis  for  stochastic 
time-variant  linear  channels  [5].  The  Input  Delay-Spread  function  h{t,  r)  maps  a 
transmitted  signal  s{t)  through  a  channel  to  the  received  signal  e{t)  as  [5,  Eq.  9]: 

=  f  h{t,t  —  T)s{T)dT. 

J  —  OO 

Fourier  transforming  h{t,T)  in  time  gives  the  Delay-Doppler- Spread  function  [5,  Eq. 
28]: 

C/(/,t)=  r  e-^^^f*h{t,r)dt. 

J-’OO 

With  each  channel  function,  Bello  associates  the  corresponding  covariance.  The  rea¬ 
son  is  that  channel  statistics  and  structure  may  be  extracted  from  the  covariance. 
For  example,  the  covariance  of  U{f,T)  under  WSSUS  assumptions  takes  the  form: 

Ruuif,  r;  /',  r')  =  E[U{r,  f)W7r)]  =  P(/,  r)5{f  -  f')5{T  -  r'), 

where  P(/,r)  =.P[|17(/,t)P]  ]s  caXLedths  scattering  function.  The  following  examples 
show  how  the  scattering  function  encodes  the  Doppler  and  delay  distributions. 

Example  1  (PMFS)  The  Phase-Modulation  Fading  Simulator  models  a  channel  as 

Hi,  '^)  =  X)  -  t„). 

Then  the  output  e{t)  is  a  sum  of  delayed  and  frequency  shifted  versions  of  the  trans¬ 
mitted  signal: 

The  Delay-Doppler  function  U (/,  r)  shows  this  “spread”  in  frequency  and  delay: 

U{f,r)  =  -^Y,an5{f  -  fn)5{T-rn). 
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If  the  (a„, /n,'r„)’s  are  IID  random  variables  with  On’s  zero  mean  and  independent 
from  (/njTn)  then  U(f,T)  is  a  generalized  random  field  with  covariance  [1,  Eq.  2.3]: 

Ruuif,  r;  /',  t')  =  (rlpf,r{f,  T)5{f  -  f')S{T  -  r')- 

Thus,  the  scattering  function  registers  the  probabifity  distribution  P/,t  (/,'’■)  of  the 
Doppler  and  delays.  The  value  of  the  PMFS  is  that  it  provides  a  simple  and  elegant 
treatment  of  the  channels  functions  and  their  covariances.  However,  there  are  non¬ 
trivial  implementation  and  ergodic  issues  with  the  PMFS  [1].  The  problems  may  be 
overcome  by  the  QMFS. 

Example  2  (QMFS)  The  Quadrature-Modulation  Fading  Simulator  modeb  a  chan¬ 
nel  as  North  &  Zeidler  [29],  Crohn  &  Bonek  [12],  Proakis  &  Salehi  [35,  pages  697- 
703]: 

N 

'^)  =  -  Tn). 

n=l 

The  output  e(t)  is  a  sum  of  delayed  and  faded  versions  of  the  transmitted  signal: 

N 

=  13  -  Tn). 

n=l 

Typically,  the  {a„(t)}’s  are  complex- valued,  JWSS  uncorrelated  Gaussian  random 
processes  with  spectrum  Pa„a„{f)  determined  by  the  propagation  environment.  If  the 
T„’s  are  also  independent  from  the  {o„(t)}’s  then 

Ruuif, t;  /',t')  =  Sif  -  /')5(r  -  r')  Pa„ar.if)PT„ir). 

n=l 

Thus,  the  scattering  fvmction  registers  the  delay  and  fading  spectrum  for  each  path. 
The  examples  show  that  the  interpretation  of  the  scattering  function  depends  on  the 
model  assumptions  or  the  propagation  environment.  For  this  reason,  an  examination 
of  the  ELB  propagation  environment  is  provided  in  Section  3.2. 


3.2  The  Surface  Scatter  Model 

For  practical  applications,  such  as  running  the  LPI  receiver  in  a  channel  emulator 
modeling  ELB  channels,  realistic  delay  and  Doppler  distributions  must  be  extracted 
from  the  channel  measurements.  This  section  applies  BeUo’s  surface  scatter  channel 
to  analyze  the  ELB  measurements.  The  resulting  framework  permits  us  to  identify 
the  QMFS  as  the  channel  model,  set  the  SNR,  the  number  of  paths,  their  delays  and 
amplitudes,  and  identify  the  fading  processes  as  Ricean.  Figure  3.1  illustrates  the 


26 


Figure  3.1:  The  surface  scatter  channel.  Prom  Bello  [1973]. 


system  setup.  The  Input  Delay-Spread  function  h{t,  r)  consists  of  three  paths  that 
modulate  the  signal  s{t).  The  signal  e{t)  at  the  antenna  is  scaled  by  the  receiver  front 
end  and  corrupted  by  additive  in-band  RF  noise  and  receiver  noise 

Figure  3.2  illustrates  the  physical  setup  [41].  A  transmitter  with  antenna  height 
hr  and  antenna  gain  Gt  is  separated  by  distance  d  from  the  receiver  with  antenna 
height  hji  and  antenna  gain  Gr.  The  direct  path  follows  free-space  propagation.  The 
coherent  reflection  also  follows  free-space  propagation  but  the  reflection  coefficient 
r  reduces  power  and  changes  phase.  According  to  Bello  [8],  anything  not  belonging 
to  these  paths  is  lumped  into  the  diffuse  scatter.  The  Input  Delay-Spread  function 
h{t,T)  is  determined  as  [24],  [9],  [26]: 

h{t,r)  -  {L5{r-  ri)  -|-  r^(r  -T2)  +g(t)5(T  -  ra)}.  (3.1) 

DirectPath  CoherentReflection  DifiuseScatter 


Figure  3.2:  Direct  path  and  coherent  reflection.  Not  shown  are  the  many  paths  that 
constitute  the  diffuse  scatter. 
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3.2.1  Direct  Path 


Basic  to  propagation  is  the  Priis  free-space  transmission  formula.  A  transmitter  with 
input  power  Py  and  antenna  gain  Gt  is  broadcasting  into  free-space  at  wavelength  A. 
At  a  distance  r,  the  receiver  intercepts  the  radiated  power  with  antenna  gain  Gr  equal 
to  the  directivity.  The  ratio  of  the  received  free-space  power  Po  to  the  transmitted 
power  Pr  is  [11],  [41,  page  68],  [10,  page  26],  [22,  page  239]: 

Go  =  ^  =  GtGr  .  (3.2) 

Additional  losses  on  the  direct  path  of  Equation  3.1  are  modeled  by  L.  The  line-of- 
sight  (LOS)  simulations  take  L  =  1.  Loo  &;  Secord  [24]  extend  L  to  a  log-normal 
process  {L{t)}  which  models  fast  fading  due  to  foliage.  North  [28]  extends  L  to 
include  multiplicative  noise  on  the  direct  and  specular  paths. 


3.2.2  Coherent  Reflection 


A  standard  factorization  for  F  is  [8],  [32,  pages  19-22],  [41,  pages  69-71]  [23,  page 
95-125]: 

r  =  DFoTZ.  (3.3) 

Divergence:  The  divergence  D  [23,  page  98]  “is  equivalent  to  a  purely  geometric 
factor  that  describes  additional  spreading  of  a  beam  of  rays  due  to  reflection  from  a 
spherical  surface”  and  has  values  between  0  and  1.  A  flat  earth  sets  D  —  1.  We  use 
the  non-linear  formula  in  Kerr  [21,  page  113]. 

Plane  Surface  Reflection  Coefficient:  Fq  denotes  the  reflection  coefficient  for  a  plane 
surface  of  the  same  material.  For  vertical  polarization,  [23,  page  99],  [10,  page  53]: 


ro  = 


esin(gg)  -  -^e  -  cosjeg)^ 
esin{6g)  +  Je-  cos{9g)^' 


where  the  grazing  angle  is 


6g  =  arctan 


hr  +  hn 
d 


and  e  is  the  complex  permittivity  of  the  terrain.  Figure  3.3  plots  Fq  for  seawater  and 
vertical  polarization  as  a  function  of  the  grazing  angle. 

Roughness  Coefficient  The  roughness  coefficient  TZ  was  developed  from  experimental 
observations  [4],  [20],  [3]  and  compares  well  with  other  work  [23,  page  122-125].  We 
use  the  modified  version  [10,  Equation  2.62]: 

^  =  e-«'/%(gV2),  (3.6) 
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Masnirude  of  Reflection  Coeffienr:  Sea  Water  20  C 


Figure  3.3:  Magnitude  and  phase  of  the  reflection  coefficient  To  for  seawater  at  20°  C 
and  vertical  polarization  as  a  function  of  the  grazing  angle  9g.  Here,  e  =  70  — j300/A. 

where  g  is  the  “electrical  roughness  parameter”  [8] 

47rA/i  sinf^o) 

g  = - (3.7) 

and  Ah  denotes  the  standard  deviation  of  the  height  fluctuation  of  the  surface.  [32, 
page  59],  Long  [23,  page  50-54]. 

Delays:  With  the  speed  of  propagation  given  by  c,  the  direct  path’s  time  delay  is 

Tx  —  c  +  {Hr  —  hxy.  (3.8) 

For  each  frequency,  the  direct  and  reflected  signals  combine  at  the  receiver.  The 
longer  path  of  the  reflected  signal  induces  a  phase  delay  A<p  given  by  [41,  page  70]: 

A0  =  27r/c(r2  -  ti).  (3.9) 


3.2.3  Diffuse  Scatter 

The  diffuse  scatter  is  modeled  by  the  zero-mean,  narrowband  complex-valued  Gaussian 
random  process  Consequently,  only  the  spectral  shape  and  variance  needs  to 

be  specified.  There  is  a  fascinating  mapping  wherein  the  real-world  propagation  man¬ 
ifests  itself  in  the  spectral  shape  [15].  Table  3.1  lists  selected  power  spectra  and  their 
associated  environments. 
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Table  3.1:  Fading  Processes 


Fading  Process 

Application 

References 

Gaussian 

Bessel 

Butterworth 

Air-to-air,  air-to-satellite 
Land  mobile  whip  antenna 
Land  mobile  to  satellite 

[8],  [9],  [16]  [26]  [40] 

[41,  page  150],  [32,  page  177] 
[24],  [25]  [25] 

Gaussian:  The  fading  process  {a(t)}  has  normalized  power  spectrum 

Paaif)  =  exp(-(//(27r/D)2) 

and  a  Gaussian  correlation  function  [8],  [9].  Here  fo  denotes  the  maximal  Doppler 
frequency  [41,  page  128]:  /p  =v/X  with  v  the  radial  speed  between  the  transmitter 
and  the  receiver  and  A  the  carrier  wavelength. 

Bessel  The  fading  process  {a(t)}  has  a  normalized  power  spectrum  [32,  page  117], 
[41,  page  150] 

p  (1  -  (//WV'"  \f\<fD 

-  I  0  elsewhere  ’ 

that  gives  a  Bessel  function  for  its  correlation.  This  fading  model  is  commonly  used  for 
land-mobile  links.  Parsons  [32,  pages  114-120]  has  a  fine  discussion  showing  that  this 
spectrum  is  the  result  of  assuming  the  multipaths  Eire  distributed  aroimd  a  vertical 
whip  antenna  uniformly  in  azimuth  but  with  zero  elevation.  Allowing  a  small  variance 
in  elevation  rolls  off  the  poles. 

Butterworth:  The  normalized  power  spectrum  is  [24] 

p  (1  +  (//*)“)-'  l/i</l> 

-  I  0  elsewhere  ’ 

and  is  used  to  model  both  multiplicative  and  the  log-normal  noise  for  land-mobile 
satellite  links  in  the  L  band. 

Power  Level  There  remains  the  problem  of  setting  the  power  level.  The  variance 
O’!  is  determined  by  the  smrface  roughness.  Beard,  Katz,  and  Spetner  [4]  used  the 
aeronautical  channel  model  to  obtain  estimates  of  cTg.  They  reported  (Tg/jroj  was 
approximately  constant  as  a  function  of  the  electrical  roughness  parameter  and  was 
bounded  between  0.2  and  0.4.  However,  later  measurements  reported  by  Beard  [3] 
showed  Cp/lPol  was  not  constant.  Although  Bello  [8]  developed  approximations  to 
this  ratio,  Ryan  ®  argues  from  conservation  laws  that  Cg  —  G/|ro|Vl  —  TV'  where  Gf 
denotes  the  “forward-scatter”  gain. 

®Frank  J.  Ryan  [1996]:  Personal  communication. 
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3.2.4  Additive  Noise 


Figure  3.4  presents  a  simplified  front  end  of  the  receiver.  The  baseband  signal  at  the 
antenna  e{t)  is  scaled  by  the  receiver’s  antenna  gain  Gr  and  cable  losses  LcaWe  at  the 
receiver  input.  The  baseband  noise  {?7(t)}  arises  from  the  RF  background  and  system 
noise  at  the  receiver  input.  Thus,  the  baseband  receiver  input  is: 

x{t)  =  GRe{t)/La>bu  + 


Figme  3.4:  Receiver  front  end. 

The  RF  background  noise  and  system  noise  is  generally  taken  to  be  zero-mean 
white  Gaussian  noise  with  fiat  spectrum  [36,  236-257],  [41,  page  310],  [9],  [40],  [13], 
[39].  By  Baseband  Assumption  B-2,  this  channel  noise  is  linearly  demodulated, 
shaped  by  the  receiver’s  matched  filter  Pfi(t),  and  then  presented  to  the  detector. 
Thus,  the  baseband  noise  {^^{t)}  is  complex- Gaussian  with  flat  spectrum  across  the 
passband  of  PR{t).  Only  the  variance  or  in-band  background  noise  power  =  cr^ 
is  needed  to  completely  specify  From  [22,  page  76]: 

P^  =  (tI  =  kT^pB,  (3.10) 

where  k  is  Boltzmann’s  constant  (1.38  x  10“^^  W  K~^  Hz“^),  B  is  the  bandwidth  of 
PR{t)  (Hz),  and  Top  denotes  the  equivalent  operating  temperature  (K).  The  literature 
on  Top  is  extensive  and  the  excellent  review  by  Pritchard  [33]  is  adapted  as  follows. 
With  respect  to  the  input  to  the  receiver,  a  simplified  model  for  a  whip  antenna  is 

Top  =  {GitTsky  + 

•^^ai»t}/Tcable  “I”  7<able(l  l/Tcable)  ”l~  T'cr* 

The  thermal  noise  sources  are  the  equivalent  temperatures  of  the  antenna  Tant,  the 
cable  Tcabie;  and  the  receiver  Tm-  The  receiver’s  antenna  gain  Gr  amplifies  the  RF 
noise  from  the  sky  T^y.  Figure  3.5  is  adapted  from  Kim  &  Muehldorf  [22,  page  74] 
and  gives  T^ky  as  a  function  of  frequency  for  RF  noise  sources. 
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Frequency:  Mm 


Figure  3.5:  Sky  temperatures  T^y  (Kelvins)  caused  by  background  noise  sources 
From  Kim  &  Muebldorf  [19951. 


3.3  ELB  Measurements 

3.3.1  Equipment  Configuration 

Figure  3.6  shows  the  configuration  of  the  mobile  transmitter,  carried  on  the  USS 
Princeton  or  the  HMMWV,  and  the  fixed-site  receiver.  A  511  pseudorandom  bit 
pattern  was  continuously  transmitted  from  the  mobile  using  BPSK  at  231.5  MHz,  40 
Watts,  999,995  bps  on  a  raised  cosine  waveform.  Over-the-air  operation  in  the  San 
Diego  area  requires  an  OA/9123  multicoupler  as  a  pre-  and  post-selection  filter  to 
prevent  transmit  jamming  of  nearby  receivers  and  to  reduce  adjacent  channel  effects 
when  received  signal  strength  is  low.  The  signal  was  sampled  using  the  AST195  at 
12  Msps  to  obtain  approximately  5.6  seconds  of  contiguous  data.  The  SABRE  GPS 
beacon  system  transmitted  once  a  minute  on  UHF  LOS  at  270.75  MHz.  The  received 
GPS  beacon  information  was  displayed  on  a  JMCIS  terminal  giving  the  location, 
bearing,  velocity,  and  past  track  information. 


Figiure  3.6:  Transmit  and  Receive  configuration. 
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3.3.2  Routes  of  the  USS  Princeton 


Figures  3.7,  3.8,  and  3.9  show  the  USS  Princeton's  location  during  the  experiment 
conducted  on  21-22  January  1998.  The  collected  channel  soundings  are  representative 
of  the  Naval  httoral  environment  for  the  Mihtary  UHF  band.  Ranges  were  extended 
because  of  a  strong  and  persistent  ducting  condition  that  existed  during  these  tests. 
A  90-nmi  range  was  obtained  by  the  SABRE  GPS  beacon  during  the  night  of  21 
January  1998  (Figure  3.8).  The  radar  range  to  the  horizon  is  approximately  35  nmi. 


"'o'  '''' 


Figure  3.8:  Position  tracks  for  the  USS  Princeton  on  21  January  1998 


Figure  3.9:  Position  tracks  for  the  USS  Princeton  on  22  January  1998 
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3.3.3  HMMWV  Routes 

Figure  3.10  maps  the  HMMWV  route  through  San  Diego  on  16  January  1998.  Prom 
1306  to  1320,  the  HMMWV  traveled  along  Catalina  Boulevard.  This  route  traverses 
a  residential  community  consisting  of  one-story  wood  houses  lining  both  sides  of 
the  two-lane  road.  Prom  1320  to  1332,  the  HMMWV  traveled  on  typical  interstate 
freeways  at  approximately  55  mph,  passing  through  downtown  San  Diego  from  1328 
to  1330.  Downtown  San  Diego  is  home  to  20  to  30  high-rises  and  numerous  multi¬ 
story  office  buildings.  Prom  1332  to  1340,  the  HMMWV  navigated  the  relatively 
small  streets  near  and  on  the  Naval  Station  San  Diego  (NAVSTASD).  NAVSTASD 
has  numerous  two-  to  four-story  concrete  buildings  close  to  the  waterfront.  Navy 
ships  line  the  piers  on  the  waterfront. 


3.3.4  An  ELB  Model 


We  close  this  section  by  comparing  the  measured  power  levels  with  the  power  levels 
predicted  by  the  model  in  Equation  3.1.  For  narrowband  signals,  it  is  standard  to 
approximate  the  effect  of  the  direct  and  reflected  paths  with  a  phase  delay.  Mathe¬ 
matically,  the  claim  is  [41]: 

-t-r5(r-T2) -|-p(t)<5(r-r3)} 

«  ^JPo{L -F  V^DTle-^^^  +  5(t)}<5(r  -  n),  (3.11) 

where  the  stationarity  of  the  {g{t)}  absorbs  the  delay.  Thus,  the  surface  scatter  model 
implies  Ricean  fading. 

How  well  does  this  model  explain  the  measmred  power  levels?  Table  3.2  hsts  the 
relevant  physical  parameters  for  the  surface  scatter  model.  The  goal  is  to  provide  an 
upper  bound  on  the  received  power.  Thus,  the  antenna  gains  are  maxiTnal  and  no 
shadowing  by  ship  or  shore  is  modeled.  Additional  assumptions  are  as  follows:  (1) 
The  direct  path  takes  L  =  1;  (2)  The  Gaussian  fading  is  small,  especially  in  bluewater, 
so  g{i)  =  0;  (3)  The  additive  noise  is  small.  Selecting  the  “Urban  Man-Made”  from 
Figure  3.5  gives  a  sky  temperatme  that  dominates  the  receiver  temperature.  The 
additive  noise  power  with  respect  to  Pr  is  101ogio(P^/Pr)  =  —149  dB.  Thus,  the 
noise  level  may  be  ignored  until  the  USS  Princeton  nears  the  horizon.  Figure  3.11 
compares  the  measured  data  to  the  free-space  and  the  surface  scatter  models.  Point 
Loma  shadows  the  transmissions  until  the  ship  rounds  the  point.  Past  this  point, 
this  model  bounds  most  of  the  data  but  does  not  capture  wide  variations  in  the 
power  levels.  Some  variation  may  be  caused  by  shadowing  of  the  transmitter  by  the 
maneuvering  ship.  Other  variations  could  be  caused  by  ducting.  Nonetheless,  that 
this  simple  model  gives  a  good  upper  bound  suggests  the  bluewater  ELB  channels  is 
amiable  to  propagation  modeling. 
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Table  3.2:  ShIp-to-Shore  parameters:  estimates  biased  high  to  get  upper  bounds  on 
received  power 


Parameter 

Symbol 

Value 

Ship:  antenna  height 

hx 

100  feet 

gain 

Gt 

OdB 

power 

Pt 

40  Watts 

Shore:  antenna  height 

hn 

355  feet 

gain 

Gr 

OdB 

Carrier  frequency 

fc 

231.5  MHz 

Receiver  Noise  Figure 

NF 

8dB 

Cable  Loss 

Lcable 

OdB 

Cable  Temperature 

Pcahle 

290  K 

Antenna  Temperature 

Psnt 

100  K 

Polarization 

vertical 

Sea  State  1 

hi/s 

0.5  feet 

Ship—to— Shore  Attenuation:  Spherical  Propagation;  Sea^State  1:  h^^\/2 


Figure  3.11:  ELB  propagation  measurements  compared  to  free-space  (dotted)  and 
the  siurface  scatter  model  (solid).  The  +’s  mark  measurements  in  San  Diego  Bay 
shadowed  by  Point  Loma.  The  additive  noise  level  is  estimated  as  10  logio( jP»j/Pt)  ^ 
-149. 
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3.4  Channel  Estimation 

Examples  1  and  2  show  that  one  must  have  a  theoretical  model  to  interpret  the  scat¬ 
tering  function.  In  practice,  the  scattering  function  must  be  estimated  from  a  finite 
number  of  discrete  samples  of  the  noisy,  distorted,  received  signal  demodulated  to 
baseband.  Consequently,  these  effects  could  confound  interpretations  of  the  scatter¬ 
ing  function.  In  our  setup,  the  digital  baseband  Delay-Spread  function  hB{m,k)  is 
smeared  in  delay  by  the  raised  cosine  filter  xjic{t)  as  [1,  Section  3.4]: 

hsim,  k)  =  xrc  *t  h{mT/Os,  kT/Os), 

where  0^=12  is  the  oversampling  rate.  The  digital  baseband  Delay-Doppler  function 
Usifik)  links  to  U{f,T)  as  [1,  Section  3.5]: 

UB{f,k)  =^XRC*r  U if,  kT/Os) 

and  the  digital  baseband  scattering  function  Pb  (/,  k)  links  to  the  scattering  function 
P{f,T)  as 

PB{f,k)=xic*rP(f,kT/0,). 

Thus,  the  channel  functions  are  smeared  by  the  raised  cosine  filter. 

Channel  estimates  are  also  corrupted  by  system  distortion,  noise,  and  demod¬ 
ulation  effects.  For  this  reason,  the  NCC  Dataset  also  contains  extensive  files  for 
calibration  and  emulated  RF  channels  [1,  Chapters  7  &  8].  Among  candidate  channel 
estimators,  an  adaptive  filter  estimate  /iAF(t,  t)  of  h{t,T)  was  found  to  be  most  robust 
against  these  effects.  It  has  a  spurious-free  dynamic  range  of  -42  dB  and  registers 
the  distortion  of  the  raised  cosine  caused  by  the  OA/9123  multicoupler.  A  DFT  of 
^Ap(^)'r)  gives  an  estimate  of  the  Delay-Doppler  function.  The  estimated  scatter¬ 
ing  function  is  obtained  by  averaging  the  magnitude  square  of  these  Delay-Doppler 
estimates. 

As  a  final  result,  QMFS  channel  parameters  are  extracted  and  tabulated  for  sub¬ 
sequent  use  in  RF  channel  emulators.  Specifically,  best  iV-path  channel  models  of  the 
form: 

N 

^n{t)S{T  -  Tn) 

n=l 

are  obtained  as  follows.  For  each  =  mT /Os,  the  2-norm  distance  between  ') 

and  5rc  •)  is  minimized  subject  to  physical  and  RF  constraints.  For  example, 

the  direct  and  reflected  paths  in  the  bluewater  ship  experiments  cannot  be  resolved 
in  time.  Thus,  only  one  path  should  be  present.  As  the  sum  of  two  passive  paths,  its 
amplitude  should  be  not  exceed  2.  More  importantly,  the  multicoupler  distorts  the 
raised  cosine.  Attempts  to  fit  a  raised  cosine  model  to  this  distorted  data  produce 
bogus  paths.  To  correct  for  this  distortion,  a  model  for  the  distorted  raised  cosine 
®Rc(*)  was  obtained  from  the  laboratory  reference  data  [1,  Chapters  7  &  8].  This 
distorted  raised  cosine  XRc(i)  model  is  used  to  fit  the  channel. 
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3.4.1  Ship-to-Shore:  Bluewater 

Data  sets  were  collected  from  the  USS  Princeton  while  departing  San  Diego  Bay  on  21 
January  1998  and  off  the  coast  of  Point  Loma  on  22  January  1998.  These  data  sets  are 
representative  of  the  Naval  littoral  environment  for  the  Military  UHF  band  at  ranges 
of  3  to  40  miles.  Figure  3.12  is  a  typical  estimate  of  the  bluewater  scattering  function. 
The  traihng  bands  in  delay  arise  from  the  raised-cosine  waveform  distorted  by  the 
OA/9123  miilticoupler.  To  see  this,  Figure  3.13  shows  the  estimated  PDP  plotted 
against  a  reference  PDP  derived  from  the  laboratory  reference  data.  The  standard 
two-path  reflection  model  predicts  a  path  delay  of  0.44  ns.  Combined  with  the  sam¬ 
pling  interval  of  83  ns,  the  two  paths  have  near-perfect  overlap.  Thus,  the  trailing 
bands  arise  from  the  distorted  raised  cosine  waveform.  With  near  perfect  overlap,  a 
single  path  model  was  fit  to  hAF(t,T)  with  the  amplitude  constraint  |ai(t)|  <2.  The 
estimated  fade  rate  fa  0.54  Hertz  is  in  concordance  with  the  maximal  Doppler 
fo  w  2.7  Hertz.  However,  this  amplitude  estimate  exhibits  a  slowly-varying  trend 
that  lowers  the  fade  rate  estimate  and  confounds  PDF  identification.  When  this  trend 
was  removed,  the  /e  wl  Hertz  and  the  detrended  amplitude  is  identified  as  Gaussian 
or  nominally  Ricean  with  K  =  (s/cr)^  =  377. 
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Figure  3.13:  Bluewater  power  delay  profile. 
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3.4.2  Ship-to-Shore:  Harbor 

In  contrast  to  this  bluewater  data,  Figure  3.14  shows  the  estimated  scattering  function 
for  the  ship  in  the  San  Diego  harbor  near  the  Coronado  bridge.  There  is  a  distinct 
second  path  delayed  by  At  «2/is  and  shifted  A/  fti-5  Hertz  from  the  direct  path. 
In  addition,  multiple  paths  extending  to  38  /xs  are  present  and  these  are  also  shifted 
by  -5  Hertz.  Figure  3.15  registers  these  delays  in  the  estimated  PDF  and  exhibits 
the  strong  second  path.  A  two-path  fit  to  h/^p{t,r)  is  also  shown  in  Figure  3.15. 
These  two  paths  are  relatively  uncorrelated  and  model  most  of  the  channel  features. 
Table  3.3  reports  the  relevant  path  parameters  as  mean  values  referenced  to  Path  1 
(/e  excepted). 

Table  3.3:  Two-Path  fit  to  Harbor;  Path  1  nominally  Ricean;  Correlation  is  0.27.  The 
"A"  denotes  the  mean  deviation  from  Path  1.  Power  normalized  to  Path  1. 


Path 

Delay 
Ar  (/xs) 

Power 

dB 

Doppler  shift 
Afo  (Hz) 

Fade  rate 
/.  (Hz) 

Fading  PDF 

1 

0.0 

0.0 

0 

1.6 

Ricean:  iC=136 

2 

-4.9 

-5 

3.8 

Ricean:  K=14.3 
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Figure  3.15:  In-harbor  power  delay  profile  and  two-path  model. 
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3.4.3  HMMWV-to-Shore:  Urban  Freeway 

Representative  of  the  HMMWV  measurements  is  the  channel  sounding  collected  on 
16  January  1998  at  13:25  while  the  HMMWV  was  driving  south  on  1-5  near  the  San 
Diego  airport  at  approximately  55  mph.  Figure  3.16  plots  the  estimated  power  delay 
profile.  The  main  lobe  is  wider  than  that  of  the  laboratory  or  ship  data  and  registers 
a  secondary  path  at  18  fis  relative  delay  implying  that  this  channel  is  dispersive 
with  respect  to  the  1-MHz  bandwidth.  Figure  3.17  plots  the  estimated  scattering 
function  and  shows  that  the  secondary  path  is  shifted  by  13  Hertz.  Moreover,  there 
are  several  paths  trailing  the  main  path.  Comparison  with  the  power  delay  profile 
and  the  delay-spread  function  shows  that  these  paths  contribute  to  the  extra  width 
of  the  main  lobe.  In  contrast  to  the  simple  two-path  model  for  the  ship,  a  multiple- 
path  model  is  required  to  capture  the  features  of  this  channel.  A  four-path  model 
was  fit  to  the  data.  Paths  1,  2,  and  3  track  the  main  lobe  while  Path  4  tracks  the 
secondary  path.  Figure  3.16  shows  that  this  model  fit  the  channel.  Table  3.4  records 
the  results.  Because  Path  2  is  slightly  stronger  than  Path  1,  the  channel  is  non¬ 
minimum  phase.  Path  4  is  relatively  weaker  and  could  be  omitted  from  simulations. 
Also,  the  corresponding  correlation  matrix  in  Equation  3.12  demonstrates  that  the 
mainlobe  paths  are  highly  correlated. 

Table  3.4:  Four-Path  fit  to  the  HMMWV;  Path  2  nominally  Ricean.  The  “A"  denotes 
the  mean  deviation  from  Path  1.  Power  normalized  to  Path  1. 


Path 

Delay 
Ar  (fis) 

Power 

dB 

Doppler  shift 
(Hz) 

Fade  rate 
f.  (Hz) 

Fading  PDF 

1 

0.00 

0.0 

0 

7 

Ricean:  K=2.d 

2 

0.98 

-fO.5 

0 

7 

Ricean:  ^=1.1 

3 

2.17 

-4.1 

1 

6 

Weibull:  K=2.0 

4 

18.16 

-20.5 

13 

4 

Ricean:  K=19 

■  1.00 

0.93 

0.71 

0.01  ■ 

0.93 

1.00 

0.89 

0.02 

0.71 

0.89 

1.00 

0.04 

0.01 

0.02 

0.04 

1.00 

(3.12) 


44 


4 


LPI  Performance  in  ELB  Multipath 


Section  2  obtained  the  LPI  performance  in  the  AWGN  channel.  Section  3  obtained 
the  fading  parameters  for  three  ELB  channels.  This  section  combines  these  results  to 
simulate  LPI  system  performance  the  ELB  channels  to  assess  multipath  effects. 

Before  we  begin  this  task,  we  need  to  extend  the  SNR  definition  to  handle  mul¬ 
tipath  and  the  LPI  receiver.  One  SNR  definition  considers  all  multipath  arrivals  as 
the  signal.  For  the  QMFS  (Example  2),  the  received  “signal”  is  therefore 

N 

=  53  -  Tn). 

n=l 

For  this  case,  the  SNR  is  the  ratio  of  crl/ag,  where  {^(t)}  is  the  additive  noise.  This 
SNR  definition  is  typically  applied  to  those  detectors  that  coherently  combine  the 
various  multipaths. 

Another  SNR  definition  considers  the  first  arrival  as  the  signal  and  the  other 
multipath  arrivals  as  interferers.  For  the  QMFS,  the  received  signal  is  then 

e{t)  =  ai{t)s{t  -  Ti) 

and  the  noise  is  g{t)  +  5^2  ~  This  SNR  definition  is  difficult  to  experi¬ 

mentally  measme  because  it  assumes  that  all  arriving  interferers  are  known  and  can 
be  measmed.  Thus,  it  is  more  appropriate  for  those  receivers  that  steer  nulls  into  the 
interferers. 

We  need  an  SNR  definition  that  directly  registers  the  multipath  effects,  permits 
immediate  and  unambiguous  scaling  to  channel  emulators,  and  also  is  amiable  to 
field  measurements.  T3rpicaUy,  noise  levels  at  the  receiver  are  acquired  by  measuring 
the  ambient  power  level  during  a  lull  in  transmissions.  With  no  fading,  the  firee- 
space  power  level  determines  the  SNR.  Therefore,  we  elected  to  use  this  SNR  in  the 
simulations.  This  SNR  also  ties  in  nicely  to  the  propagation  factor  for  real-world  mea¬ 
surements.  For  example,  if  there  was  a  credible  RF  model  for  propagation,  then  firom 
the  measured  noise  levels,  the  RF  mission  planner  could  select  a  transmission  power 
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Pt  to  set  an  operational  SNR.  Computationally,  this  means  that  the  AWGN  simu¬ 
lations  of  Section  2  are  replicated  with  the  following  difference.  Instead  of  receiving 
y{t)  =  s{t)  +  g{t),  we  receive 


N 

y{^)  =  Ylan{t)s{t-Tn)  +  g{t). 

n=l 

Path  1  is  the  first  arrival  with  power  Pi  =  £^[|ai(t)p].  The  SNR  is  set  as 

SNR  = 

The  relative  power  of  the  remaining  paths  are  scaled  to  Path  1.  No  path  correlations 
are  assumed  in  these  simulations. 

This  section  runs  the  LPI-1  and  LPI-2  simulations  using  the  AWGN  level  of 
Section  2.  The  Bluewater,  Harbor,  and  HMMWV  multipath  simulations  of  Section  3 
are  scaled  by  our  SNR  definition.  Thus,  any  deviations  firom  the  AWGN  simulations 
are  caused  by  the  multipath.  Finally,  all  simulations  assume  a  chip  duration  of  rc=10 
ys. 

4.1  LPI-1  ELB  Multipath  Performance 

For  the  LPI-1  simulations,  the  fading  processes  are  Ricean  with  Gaussian  spectral 
shape.  Referring  to  Section  3.2.3,  this  means  that 

an{t)  =  Pn  X  {Kn  +  Pn(t))) 

where  {pn(t)}  is  a  zero-mean,  unit-variance,  complex-valued,  WSS,  Gaussian  random 
process.  Its  spectrum  is  proportional  to  exp((//(27r/D))^)  with  fo  set  to  the  fade 
rate  reported  for  the  Bluewater,  Harbor,  and  HMMWV  multipath  simulations.  The 
Kn  is  also  set  to  the  Ricean  K  factor  for  these  simulations.  The  power  level  is 
Pn  =  £^[|an(t)p]  and  rigged  so  that  the  relative  power  levels  scale  to  the  dB  power 
levels  reported  for  the  Bluewater,  Harbor,  and  HMMWV  multipath  simulations. 

4.1.1  Performance  for  [Np,  iVJ  =  [1024, 0, 2048] 

Figures  4.1,  4.2,  and  4.3  compare  the  performance  of  LPI-1  in  AWGN  and  the  mul¬ 
tipath  simulations.  The  LPI-1  setup  is  the  low-dimension  simulation  used  in  Sec¬ 
tion  2.2.2.  The  PN  sequence  has  length  N^.  =  2048  and  is  partitioned  to  give  L  =  2 
symbol  vectors  of  length  Np  =  1024.  With  the  chip  duration  set  to  rc=10  gs,  then 
the  symbol  vector  has  dmration  10  ms.  With  all  fade  rates  less  than  5  Hertz,  there  is 
little  fading  variation  on  each  symbol  vector. 
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The  fading  is  small  in  the  Bluewater  multipath.  Consequently,  Figure  4.1  shows 
little  difference  between  the  AWGN  and  the  multipath.  The  Harbor  multipath  has  two 
paths.  The  second  path  is  delayed  by  At  =  2  fjs  and  5  dB  down  from  the  first  path. 
Figure  4.2  shows  that  the  Harbor  multipath  reduced  the  probability  of  symbol  error 
by  10  dB.  In  contrast,  Figure  4.3  shows  HMMWV  multipath  significantly  increased 
the  probability  of  symbol  error.  Of  the  four  HMMWV  paths.  Path  1  and  Path  2 
both  have  significant  fading,  similar  power  levels,  and  a  delay  of  only  At  =  0.98  fis. 
The  large  K  factors  of  the  Ricean  paths  in  the  Harbor  multipath  result  in  a  degree 
of  temporal  diversity  that  enhances  detector  performance.  The  smaller  K  factors  in 
the  HMMWV  multipath  imply  deeper  fades,  resulting  in  a  symbol  error  floor. 


Figure  4.1:  LPI-1  [Np,  N^,  =  [1024,0,2048]  performance  in  Bluewater  multipath. 
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PN(ll^l);  Partial  Sequence  Length  L=2 


Figure  4.2:  LPI-1  [iVp,  N^,  iVu]  =  [1024, 0, 2048]  performance  in  Harbor  multipath. 


Figure  4.3:  LPI-1  [Np,  N^,Nu]  =  [1024, 0, 2048]  performance  in  HMMWV  multipath. 
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4.1.2  Performance  for  [A^p,  ^u]  =  [1024, 0, 2^®] 

Figures  4.4,  4.5,  and  4.6  compare  the  performance  of  LPI-1  in  AWGN  and  the  mul¬ 
tipath  simulations.  This  LPI-1  setup  follows  the  setup  of  Section  2.2.3.  The  PN 
sequence  has  length  Nu  =  2^®  and  is  partitioned  to  give  L  =  64  symbol  vectors  of 
length  Np  =  1024.  The  simulations  replicate  the  results  just  obtained  in  Section  4.1.2: 
no  difference  between  Bluewater  multipath  and  AWGN;  Harbor  multipath  improves 
performance;  HMMWV  multipath  degrades  performance.  Thus,  the  number  of  sym¬ 
bol  vectors  do  not  affect  multipath  performance. 


Figure  4.4:  LPI-1  [iV'p,i\rA,iV’„]  =  [1024,0,2^®]  performance  in  Bluewater  multipath. 
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SNR  (dB):  Nuiober  of  suDulations=  100000;  G:rautan  spednn 


Figure  4.5:  LPI-1  [iVpjiVAjiV,/]  =  [1024,0,2^®]  performance  in  Harbor  multipath 


Figure  4.6:  LPI-1  [Np,N^,Nu\  =  [1024,0,2^®]  performance  in  HMMWV  multipath 


4.2  LPI-2  ELB  Multipath  Performance 

Section  4.1  showed  that  the  number  of  symbol  vectors  do  not  affect  multipath  perfor¬ 
mance.  For  this  reason,  the  LPI-2  simulations  were  undertaken  to  see  how  the  length 
of  the  symbol  vectors  affected  multipath  performance.  For  these  LPI-2  simulations, 
the  fading  processes  are  Ricean  with  Bessel  spectral  shape  described  in  Section  3.2.3. 
Simulations  showed  that  for  these  relatively  slow  fade  rates,  changing  the  spectral 
shape  between  Gaussian  and  Bessel  did  not  significantly  affect  the  probabihty  of 
symbol  error. 

The  LPI-2  setup  is  the  same  as  in  AWGN  simulations  in  Section  2.3.  The  PN 
sequence  has  length  =  2^®.  The  CCSK  symbol  vectors  have  lengths  Np=204l8, 
8192, 16384,  and  65536.  With  the  chip  duration  set  to  Tc=lQ  /is,  the  symbol  vectors 
have  durations  20,  82, 164,  and  665  ms.  Although  all  fade  rates  are  less  than  5  Hertz, 
these  longer  symbol  vectors  can  have  fading  variations. 

Figures  4.7,  4.8,  4.9,  and  4.10  plot  the  LPI-2  performance  for  symbol  vectors  of 
length  iVp=2048,  8192,  16384,  65536.  Each  plots  shows  the  AWGN,  Bluewater,  Har¬ 
bor,  and  HMMWV  simulations.  As  expected,  the  Bluewater  multipath  shows  little 
difference  from  the  AWGN.  For  the  Heirbor  multipath,  LPI-2  with  the  short  partial 
lengths  behaves  like  LPI-1  with  better  than  AWGN  performance.  However,  as  the 
partial  length  increases  then  LPI-2  performance  degrades  to  AWGN  performance.  For 
HMMWV  multipath,  LPI-2  with  short  partial  lengths  shows  worse  than  AWGN  per¬ 
formance.  However,  increasing  the  partial  length  slowly  improves  LPI-2  performance. 
Past  a  threshold,  there  is  a  sudden  improvement  to  better  than  AWGN  performance. 

So,  for  LPI-2  to  get  better  than  AWGN  performance,  we  should  select  short  partial 
lengths  in  the  Harbor  and  a  long  partial  length  in  the  HMMWV.  This  is  not  imme¬ 
diately  obvious  why  this  is  so.  This  contradictory  performance  backs  the  assertion 
in  Section  1.1  that  the  LPI  system  also  include  an  operational  context.  Moreover, 
the  partial  length  is  only  one  “knob”  of  the  [Np,N/^,Nu]  parameters.  A  non-trivial 
question  is  to  determine  optimal  LPI  performance  in  multipath  by  adjusting  all  three 
parameters.  Addressing  this  question  constitutes  the  bulk  of  the  Phase  4  recommen¬ 
dations  listed  in  Section  1.1. 
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Figure  4.7:  LPI-2  [Np,N^,  Nu]  =  [2048,2047,2^®]  performance  in  multipath 


Figure  4.8:  LPI-2  [ATp,  iV^,  iV«]  =  [8192,8191,2^®]  performance  in  multipath. 
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